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Abstract 

The problem of evaluation of Lyapunov exponent in queueing net- 
work analysis is considered based on models and methods of idem- 
potent algebra. General existence conditions for Lyapunov exponent 
to exist in generalized linear stochastic dynamic systems are given, 
and examples of evaluation of the exponent for systems with matrices 
of particular types are presented. A method which allow one to get 
the exponent is proposed based on some appropriate decomposition of 
the system matrix. A general approach to modeling of a wide class 
of queueing networks is taken to provide for models in the form of 
stochastic dynamic systems. It is shown how to find the mean service 
cycle time for the networks through the evaluation of Lyapunov expo- 
nent for their associated dynamic systems. As an illustration, the mean 
service time is evaluated for some systems including open and closed 
tandem queues with finite and infinite buffers, fork-join networks, and 
systems with round-robin routing. 



1 Introduction 

The evolution of actual systems encountered in economics, management, en- 
gineering, and other areas can frequently be represented through stochastic 
dynamic models of the form 

x(k) = A(k)x(k - 1), 

where A(k) is a random state transition matrix, x(k) is a state vector, and 
matrix-vector multiplication is thought of as defined in terms of a semiring 
with the operations of taking maximum and addition [IJ [2j [31 HI [5] . 

In many cases, the analysis of a system involves evaluation of asymptotic 
growth rate of the system state vector x{k) , which is normally referred to as 
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Lyapunov exponent [21 [6]. In the semiring with the operations of maximum 
and addition, the mean growth rate is defined in terms of the semiring as 

A= lim HsbCAOH 1 /*. 

k— too 

Note that for queueing networks, the value of A can be considered as 
the mean service cycle time in a system, while its inverse can be thought of 
as system throughput. 

Existence conditions for the above limit can normally be established 
based on the ergodic theorem in [7j. Specifically, examples of the conditions 
can be found in [H El [8]). Note however, that these conditions frequently 
require that the system state transition matrix is of particular form or type 
(matrix with finite entries, irreducible matrix). 

For systems with fixed nonrandom matrices, the value of A can easily 
be found based on results in 0[TO] (see also OUT]. However, for stochastic 
systems with random state transition matrices, evaluation of A normally 
appears to be a difficult problem. Existing results [H Q21 El Q31 [H| include 
the exact solution to the problem for second-order systems determined by 
matrices with entries that are independent and have either normal or ex- 
ponential probability distributions. It is shown in |15l [TBI [T71 [T8] how the 
value of A can be found for systems with particular matrices of arbitrary 
size, including triangular matrices, similarity matrices, and matrices of rank 
1. 

In the current paper we consider the problem of evaluation of Lyapunov 
exponent, which arises in the analysis of a wide class of queueing networks. 
We start with an overview of some concepts and results in idempotent alge- 
bra that underlie subsequent parts, and introduce related notations. 

We consider a general model of stochastic dynamical system governed by 
the vector equation with the random matrices A(k) which are assumed to 
be independent and identically distributed for all k = 1,2,... New general 
conditions for the limit that determines the Lyapunov exponent A to exist 
are given. 

We present results of evaluating A for systems with matrices of particular 
types, including diagonal and triangular matrices, similarity matrices, and 
matrices of rank 1. A new approach to evaluation of A is then proposed 
based on a decomposition of A{k) . The approach allows one to reduce the 
problem with the original matrix A(k) to that with a matrix of particular 
type and known solution. 

Furthermore, we discuss application of the above results to the analysis 
of a class of queueing systems described in [191 1201 IZU [22] . Some particular 
systems including open and closed tandem queues with finite and infinite 
buffers, a fork-join network, and a system with round-robin routing are 
examined. For these queueing systems, we obtain the mean service cycle 
time through evaluation of Lyapunov exponent in their associated dynamical 
systems. 
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2 Preliminary results 



2.1 Idempotent algebra 

Consider a set X with two operations © and <8> referred to as addition and 
multiplication, and neutral elements and 1 called zero and identity. We 
suppose that (X,©, <g>) is a commutative semiring with idempotent addi- 
tion and invertible multiplication. Such a semiring is normally called an 
idempotent semifield. 

Let X + = X \ {0}. Then each x G X + is assumed to have the inverse 

. For any x, y G X + , the power x y is defined in an ordinary way. Note 
that in what follows, the sign <g> will be omitted as is usual in conventional 
algebra. The notation of power is thought of as defined in terms of idempo- 
tent algebra. However, for the sake of simplicity, we use ordinary arithmetic 
operations in the expressions that represent exponents. 

Since the addition is idempotent, one can define a linear order < on 
X according to the rule: x < y if and only if x © y = y . Below the 
relation symbols < will be understood only in the sense of this linear order. 
According to the order, for any x G X, we have x > 0. 

Examples of idempotent semirings (semifields) include 

Kmax,+ = (RU {-oo}, max, +), (R min ,+ = (R U {+oo}, min, +), 
Kmax.x = (R+ U {0},max, x), R min>x = (R+ U {+oo},min, x), 

where R is the set of real numbers, R + = {x G R\x > 0}. 

Consider the semiring R maXj+ . It is easy to see that = — oo, 1 = 0. 
For each x G R, there exists its inverse x~ l , which is equal to —x in the 
ordinary arithmetics. For any x,y G R, the power x y coincides with the 
arithmetic product xy . The linear order has conventional meaning. 

For the semiring Rmin,x > we have — ©oo, 1 — 1. The inverse element 
and the power have ordinary meaning. The relation < determines an order 
that is the reverse of the natural linear order on R + . 

It is easy to verify that the semirings R max ,+ , K m in,+ , K maXjX , and R min ,x 
are all isomorphic. 

2.2 Matrix algebra 

For any matrices A,B G X mx " , C G X nx ' , and a scalar x G X , the matrices 
A © B , BC , and x A are defined in a usual way with the formulas 

n 

{A(BB}ij = {A} tj e{B} ij7 {BC}ij = Q){B} lk {C} kj , {xA} tj = x{A} iy 

k=i 

The matrix with all its entries equal to is referred to as zero matrix 
and denoted by 0. A matrix is called regular if every row of the matrix 
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has at least one nonzero entry. Below we assume all nonzero matrix to be 
regular. 

For each matrix A = (aij) € X mxn , its norm is defined as 

n m 

W A W = (D ® ai i- 

i=l j=l 

For any matrices A and B of appropriate size, and a scalar x > 1, it 
holds 

||Ae5|| = p|| e||B||, \\AB\\ < ||A||||J3||, \\xA\\ = x\\A\\. 

A square matrix is called diagonal if all its off-diagonal entries are zero, 
and triangular if its entries either above or below the diagonal are zero. The 
matrix / = diag(l, . . . , 1) is referred to as identity matrix. The integer 
nonnegative powers of a matrix A are defined by the formulas A = I, 

A k+i = A k A i for all k,l = 0,l,... 

A matrix A E x nxn is a matrix of a similarity operator (similarity 
matrix) if there exists a constant a/C called the coefficient of similarity, 
such that for any x € X n , it holds 

||Aa;|| = q 1 1 £c 1 1 . 

A vector y £ X n is said to be linearly dependent of x±,... ,x m G X n 
if y = aixi © • • • © a 

m-Em for some scalars oi, . . . ,o m S X. The rank of a 
matrix A is its maximum number of linearly independent rows (columns). 
A matrix A has the rank 1 if and only if A = xy T , where x and y are 
some nonzero vectors. 

The sum of diagonal entries of a matrix A = (aij) E X nxn is called its 
trace, and denoted by 

n 

tr A = cijj. 
i=l 

A number A is called an eigenvalue of the matrix A if there exists a 
vector x ^ such that Ax = Aa; . 

The maximum (in the sense of the linear order on X ) eigenvalue is called 
the spectral radius of A and evaluated as 

n 

p(A) = ®tr 1 / m (A m ). 

m=l 

The next result has been obtained in 0[TU] (see also [51 ITT]). 
Theorem 1 (Romanovskii) . For any matrix A G X nxn there exist limits 
lim || A* H 1 /* = p(A), lim tr 1/k (A k ) = p(A). 

fc— >oo k— ¥oo 
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2.3 Properties of the expected value 



Consider some inequalities which relate evaluation of the expected value of 
random variables to operations involved in idempotent semirings. To be 
definite, assume that the semiring of interest is R m ax,+ • 

We suppose that all random variables under considerations are defined 
on a common probability space and have finite mean values. 

It is easy to verify that for any random variables £ and r] it holds 

E(£ © rj) > E £ © E V, E£v = E £E?7- 

Let A be a random matrix. We use the symbol EA to denote the matrix 
obtained from A by replacing all its entries with their expected values, 
provided that EO = . 

For any random matrices A and B of appropriate size, it holds 

E(A © B) > EA © EB, EAB > EAEB, E\\A\\ > \\EA\\, EtrA > tr(EA). 

3 Stochastic dynamical systems 

Let A{k) € x nxn be a random state transition matrix, x{k) £ X n be a state 
vector, k = 1, 2, . . . Consider a dynamical system described by the equation 

x(k) = A T (k)x(k - 1). 

Note that the representation in the form with the transpose is intended 
to simplify further formulae. Clearly, the use of the matrix A T (k) does not 
actually change the general form of the equation which is normally written 
as x(k) = A(k)x(k - 1). 

In particular, the above representation allows products of matrices 
A(l), . . . , A(k) to be examined in the natural order. With the notation 

A k = A(l)..-A(k), 

iterating the dynamic equation gives x{k) = A£x(0). 

Suppose that the sequence {A(k)\k > 1} consists of independent and 
identically distributed random matrices, E||Ai|| is finite. However, we do 
not require that for each k = 1,2, ... , the entries of A(k) are independent. 

Note that most of the statements below remain valid if the matrices 
A(k) form a stationary sequence rather than a sequence of independent and 
identically distributed matrices. 

3.1 Lyapunov exponent 

In many cases, the analysis of a system involves evaluation of asymptotic 
growth rate of the system state vector x(k) , which is frequently referred 
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to as Lyapunov exponent. In the semirings IR m ax,+ and K m in.+ ) the mean 
growth rate is defined as 

A = lim \\x(k)\\ 1 / k (1) 

fc— too 

provided that the limit on the right side exists. 
In R maXj x and K m i n ,x , we have 

A= lim log\\x(k)\\ 1/k . 

fc— >oo 

In what follows, we will concentrate on the semiring R max ,+ • Taking into 
account isomorphisms between R max ,+ and the semirings R m i ni +, K m ax,x 
and R m in,x , all results below can easily be extended to these semirings. 

Suppose that the entries of the initial vector £c(0) are finite with prob- 
ability one (w. p. 1). Then the mean growth rate A can be defined as 

A= lim p fc || 1/fc . (2) 

fc — >OD 



3.2 Existence conditions 

General conditions for limit ([2]) to exist can be established based on the clas- 
sical result in [7j, which is presented below in terms of ordinary arithmetic 
operations. 

Theorem 2 (Kingman). Let {CimV < m l ^ e a family of random variables 
which satisfy the following properties: 

1. (im< Cik + Cfcm (subadditivity); 

2. The joint distributions are the same for both families {CimV < Tn} and 
{Cz+i, m +i|^ < m } (stationarity); 

3. For all n = 1,2,..., there exists E£ofc > — ck for some positive con- 
stant c (boundedness) . 

Then there exists a constant X, such that it holds 

lim C ok /k = A w. p. 1, lim E( Qk /k = A. 

k—>oo k— >oo 

For the semiring (R max ,+ > the existence conditions can be formulated as 
follows. 

Theorem 3. Let {A(k)\k > 1} be a stationary sequence of random matri- 
ces, E||vli|| < oo and p(E[Ai]) > — oo. Then there exists a finite number A 
such that 

lim \\A k \\ 1/k = A w. p. 1, lim E||A fc || 1/fc = A. 

fc— >oo fc— too 
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Proof. Consider the family {(i m \l < m} with Q m = \\A(l + l)A(l + 
2) • • • -A(m)|| , and note that Cok = \\Ak\\ ■ 

Let us verify that the conditions of Theorem [2] are fulfilled. 

Subadditivity follows from the inequality 

\\A(l + 1) • • • A(m)\\ < \\A(l + 1) • • • A(k)\\\\A(k + !)■■■ A(m)\\. 

It is clear that the family is stationary. From the condition E||Ai|| < oo 
it follows that E||Afe|| < oo. Furthermore, since p(EAi) > = — oo and 
E\\A k \\ = E\\A(1) ■ ■ ■ A(k)\\ > ||(EAi) fc || > p k (EA 1 ), the family is bounded. 

By applying Theorem [21 we arrive at the desired result. □ 

4 Evaluation of Lyapunov exponent 

Now we give examples of evaluation of Lyapunov exponent for systems with 
matrices of particular types [151 US EZ1 EE] • Note that for each system under 
consideration, the conditions of Theorem [3] are assumed to be fulfilled. 

4.1 Systems with diagonal matrix 

Suppose that A(k) = diag(<ii(fc), . . . , d n (k)) is a diagonal matrix, k = 
1,2,... It is easy to verify that 

A = tr(EAi). 

4.2 Systems with similarity matrix 

Let A(k) be the matrix of a similarity operator, k = 1, 2, . . . Then it holds 

A = E\\Ax\\. 

4.3 Systems with matrix of rank 1 

Suppose that for all k = 1,2,..., there exist vectors u(k) and v(k) such 
that A{k) = u(k)v T (k). Then it holds 

A = E[/(1W2)]. 

4.4 Systems with triangular matrix 

The following result was obtained in [18J based on classical bounds for the 
mean values of sums of independent random variables combined with new 
inequalities for products of triangular matrices in idempotent algebra. 

Theorem 4. If the matrix A(k) is triangular k = 1,2,..., then it holds 

A = tr(EAi). 
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4.5 A matrix decomposition method 

If the state transition matrix of a system falls into the one of the particular 
types examined above, evaluation of Lyapunov exponent presents no special 
problems. In the case that the matrix has a different type, one can try to 
implement the following approach. 

Let us assume that there exists a decomposition of the matrix A(k) in 
the form 

A(k) = B(k)C(k), 

where B(k) and C(k) are some independent matrices. If a solution to the 
problem of evaluation of Lyapunov exponent for the system with the matrix 

A'(k) = C(k)B(k + 1) 

is known, it can normally be taken as the solution for the system with the 
initial matrix A(k) . Otherwise, one can continue with decomposition of 
A'(k). 

Note that for any matrix A(k) of rank 1 , the above decomposition exists 
and takes the form A(k) = u(k)v T (k). 



4.6 A system with a matrix of incomplete rank 

Consider a matrix A of order n. Suppose that rank A < n. Then there 
exits a skeleton decomposition of the matrix 

A = BC. 

The above decomposition will be referred to as backward triangular if the 
matrix CB is triangular. 

Theorem 5. If the matrix A{k) allows for the backward triangular skeleton 
decomposition A(k) = B{k)C{k) with independent factors B(k) and C(k), 
then 

A = tr(E[C(l)B(l)]). 

Proof. Consider the matrix 

k ( k-i \ 

A k = II B (j)C(j) = B(l) I] C(j)B(j + 1) C(k). 
3=1 \i=i / 

Furthermore, we have 

fc-i 



E\\A k \\ > EtrA k > E 



tr C(k)B(l) ft C(j)B(j + 1) 



>tr(E[C(l)£(l)] fe ) 
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By applying Theorem [TJ we get 



lim EPfcH 1 ^ > lim ( tr(E[(7(l)B(l)] fc ) 



1/fc 



p(E[C(l)B(l)]). 

Since C(l)fi(l) is a triangular matrix, we have the inequality 

A>tr(E[C(l)B(l)]). 

It remains to verify that the opposite inequality is also valid. First we 
write 

fc-i 

||4b||< l[C(j)B(j + l) \\B(l)\\\\C(k)\\. 

3=1 

The matrices C(j)B(j + 1) are triangular and independent for all j = 
1, . . . , k — 1 . With Theorem HJ we get 

A<tr(E[C(l)B(l)]). □ 



5 Algebraic models of queueing networks 

Queueing networks with fork-join operations present a quite general class of 
dynamical systems that can be described in terms of the semiring (R max ,+ 
by the equation [EU |2JB EH |22J E3] 



su(fc) = A(k)x(k - 1). 



(3) 



The fork and join operations allow customers (jobs, tasks) to be split into 
parts, and to be merged into one, when they circulate through the system. 
The fork-join formalism proves to be useful in the description of dynamical 
processes in a variety of actual systems, including production processes in 
manufacturing, transmission of messages in communication networks, and 
parallel data processing in multi-processor computer systems. As an illustra- 
tion of the fork and join operations, one can consider respectively splitting a 
message into packets in a communication network, each intended for trans- 
mitting via separate paths, and merging the packets at a destination node 
of the network to restore the message 



5.1 Fork-join queueing networks 

Consider a network with n single-server nodes and customers of a single 
class. The network topology is described by an oriented graph Q = (V, E) , 
where V = {1, . . . ,n} is a set of nodes, and E = C V x V is a set 

of arcs that determine the transition routes of customers. 

For every node i G V, we define the sets P(i) = € E} and 

= ^ E} ■ The nodes i with P(i) = are assumed to be 
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source nodes that represent infinite external arrival streams of customers. 
If S(i) = 0, the node i is considered as an output node intended to release 
customers from the network. 

Each node i includes a server and a buffer which together present a 
single-server queue operating under the first-come, first-served (FCFS) dis- 
cipline. At the initial time, all servers are free of customers. The buffer at 
each source node has infinite number of customers, whereas the buffer at 
any other node i has q customers, < q < oo . 

We suppose that in the network, the usual service procedure is combined 
with additional join and fork operations [23] which may be performed in a 
node respectively before and after service of a customer. The join operation 
is actually thought to cause each customer which comes into node i, not to 
enter the buffer at the server but to wait until at least one customer from 
every node j £ P(i) arrives. As soon as these customers arrive, they, taken 
one from each preceding node, are united to be treated as being one customer 
which then enters the buffer to become a new member of the queue. 

The fork operation at node i is initiated every time the service of a cus- 
tomer is completed; it consists in giving rise to several new customers instead 
of the former one. As many new customers appear in node % as there are 
succeeding nodes included in the set S(i) . These customers simultaneously 
depart the node, each being passed to separate node j € S(i) . We assume 
that the execution of fork-join operations when appropriate customers are 
available, as well as the transition of customers within and between nodes 
require no time. 

5.2 Dynamical equation 

Let Tik be the kth service time, and Xi(k) be the kth departure epoch at 
node i = 1, ... ,71. 

We assume that for all i = 1, . . . , n, the sequence {r^, k > 1} consists of 
independent and identically distributed nonnegative random variables with 
a finite mean value. Furthermore, for each k , the random variables and 
Tjk are independent for all i ^ j . 

With the condition that the network starts operating at time zero, we 
put Xj(0) = = 1 and Xi(k) = — oo = for all k < 0. We introduce 
notation 





( 


xi(k) \ 




( Tl k 





x(k) = 














{ 


x n (k) J 






Tnk 



Let us define M = max{cj|cj < oo, i = 1, ...,n}. For each m = 
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0,1, . . . , M , we introduce the matrix G m = (<7^?) with entries 

m (l, if i G P(i) and m = Cj, 
%3 |0, otherwise. 

It is easy to see that the matrix G m can be considered as an adjacency 
matrix of the partial graph Q m = (V,E m ), where E m = G P(j),Cj = 

m} . 

Now we can formulate the following result [201 122] [23] . 

Lemma 6. Suppose that the graph Qq associated with the matrix Gq is 
acyclic, and r is the length of its longest path. Then the dynamics of the 
network is described in the semiring R max ,+ by the equation 

M 

x(k) = Q)A m (k)x(k-m), (4) 

m=l 

where 

A 1 {k) = (I@T k Gl) r T k {I®G T l ), 
A m (k) = (I@T k Gl) r T k Gl l , m = 2,...,M. 

5.3 Networks with finite buffers 

Suppose now that the buffers at servers in the network may have limited 
capacity. In such a network, servers may be blocked according to some 
blocking mechanism [241 [20] . Below we consider networks operating under 
the manufacturing and communication blocking rules which are commonly 
encountered in practice. 

Manufacturing blocking at node i implies that a customer cannot release 
the server if there is at least one succeeding node j G S(i) without empty 
space in its buffer. The communication blocking rule requires the server in 
node i not to initiate service of a customer until there is an empty space in 
the buffer in each node j G S{i) . 

Suppose that the buffer at node i has capacity b{ , < bi < 00 . Clearly, 
for all i = 1, . . . , n , we have bi > Cj . 

Let us define M\ = max{cj|cj < 00, i = 1, . . . ,n} and M2 = max{6j|6j < 
00, i = 1, . . . , n} + 1 . For each m = 1, . . . , M2 , we introduce the matrix 
H m = (h^-) with its entries 

, m= fl, if j G S(i) and m = bj + 1, 
^ [O, otherwise. 

Consider M = max{Mi,M 2 }. If M 2 > Mi, then put G m = for all 
m = M l + l,M l + 2,...,M 2 . 

We have the following result [20] [22] [23] . 
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Lemma 7. Suppose that the graph Qq associated with the matrix Gq is 
acyclic, and r is the length of its longest path. Then the dynamics of the 
network is described in the semiring K m ax,+ by the equation Q with the 
matrices defined under manufacturing blocking rule as 

A l (k) = (I © T k G^) r (T k {I © Gj) © ill), 
A m (k) = (I®T k Gl) r {T k Gl l @H m ), m = 2,...,M; 

and under communication blocking rule as 

Ai(k) = {I®T k Gl) r T k {I@Gl @H X ), 
A m (k) = (I © T k G%) r T k (Gl ®H m ), m = 2, . . . , M. 

6 The mean service cycle time in queueing sys- 
tems 

Consider a queueing system and suppose that its dynamics is described by 
equation Q. The evolution of the system can be represented as a sequence of 
service cycles. The first cycle starts at the initial time, and it is terminated as 
soon as all the servers in the network complete their first service, the second 
cycle is terminated as soon as the servers complete their second service, and 
so on. Clearly, the completion time of the kth cycle can be represented as 
||a?(fc)|| . With the condition £c(0) = 0, we have 

||x(A0|| = pfc||, A k = A(k)---A(l). 

In the queueing systems, the mean growth rate of the state vector (Lya- 
punov exponent) can be regarded as the mean service cycle time. The recip- 
rocal of the mean growth rate can be considered as the system throughput. 

Below we show how for some queueing networks, the mean service cycle 
time is evaluated based on examination of their related system matrices 
A(k). 



6.1 Open and closed tandem queues 



Tandem queueing systems present networks with the simplest topology de- 
termined by graphs which include only the nodes with no more than one 
incoming and outgoing arcs. Consider an open system with n nodes and 
infinite buffers, depicted in Figure [TJ 

Since M = , we put G\ = and then get equation ([3|) with the matrix 



A(k) = (/ © T k G\ 



T\n-1 



n 



Tlk 
T\kT2k 



T~2k 



\ T lk ■ ■ ■ T nk T 2k 



\ 



Tnk ) 
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1 



2 



Cl = OO C2 = 



Figure 1: Open tandem queues with infinite buffers. 



The state transition matrix is triangular. The application of Theorem 2] 
yields 

A = tr(ETi) = max(Ern, . . . , Er„i). 
Consider a closed tandem queueing system presented in Figure [2j 



1 2 

OO — OO 



71 



Cl 



^2 



Figure 2: A closed tandem queueing system with infinite buffers. 

Suppose q = 1 for all i = 1, . . . ,n. Then we have M = 1. The state 
transition matrix takes the form 



A{k)=T h {I®Gl) 



( r lk 

T~2k T~2k 



V 



With n = 2, the matrix is a similarity matrix with the factor 

|A(fe)|| = Tifc © t 2fc = trTfe. Therefore, we have 

A = EtrTi = E max(rn, T2i). 



6.2 Tandem queues with finite buffers and blocking 

Consider an open tandem queueing system which may have finite buffers 
(Figure [3|). 

Suppose that the system operates under the manufacturing blocking rule. 
Let n = 3 , b\ = 62 = 00 , and 63 = . Then we have M = 1 . Evaluation of 
the system matrix gives 



A(k) = (/er fc G^) 2 (r fc ©^i 



Tlk v v 

TlkT2k T 2 k 1 
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1 2 n 

□O — 00 — * 00 — * 

Cl = OO c 2 = c„ = 

bi = oo 6 2 6 n 

Figure 3: Open tandem queues with finite buffers. 

Taking into account collinearity of the last two rows, we get the skeleton 
decomposition 

Tlk 

A(k) = ( Ti k T 2k r 2k 1 

TlkT~2kT~3k T2kT3k T 3k 

o i ) ( Tik - : i = B(k)c{k). 

1 1 r lk r 2k r 2k 1 1 



T 3k 



Consider the matrix 



C(k)B(k + l)= Tlk ° 

V T ikT 2k T 2k © r 3ifc+ i 

By application of Theorem [5l we conclude that 

A = tr(E[C(l)B(l)]) = max(Ern, Emax(r 2 i,r 3 i)). 

Assume that the system in Figure [3] follows the communication blocking 
rule. With n = 3, b\ = oo, and b 2 = b 3 = 0, we have M = 1, Gi = 0, and 
Hi = Go • The state transition matrix is represented as 

/ rik nk <D 

A(k) = (I®T k GlfT k {I®G G ) = r lfe r 2fc r lfc r 2fe r 2fc 

\ TikT 2k T 3 k TlkT 2 kT 3 k T 2k T 3k 

Furthermore, we get the decomposition 

1 

•is/.) - I r 2fe r 2fe I ( T i* T i fe ;; ) Da-)cui 

T2kT3k T 2 kT 3 k 

Consider the matrix product 



C(k)B(k + l) 



TlkT2,k+l T lkT2,k+l 
T2,k+lT 3 ,k+l 72,fc+l T 3,fc+l 

Tlk 
T3,k+1 



T2,k+i T 2 ,k+i ) = u(k)v T {k). 
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Since the product is actually a matrix of rank 1 , we finally get 
A = E[v T (l)u(2)] = Er 2 i + Emax(ni,T 3 i). 

6.3 A fork- join network 

Now we turn to network models which involve fork and join operations. 
Consider a fork-join network with n = 5 nodes and infinite buffers, depicted 
in Figure HI [22]. 




c 4 = 



Figure 4: A fork-join queueing network with infinite buffers. 
Since r = 2, we arrive at equation ([3]) with the matrix 
A(k) = (I@T k GlfT k {I@G T 1 ) 



( T lk 



\ 





T2kT3k T 2 kT~3k «J 

T 3k T 3k 

TXkTlkTAk T2 k T 3k T ik T2 k T 3k T ik T± k 

T 5k T 5k T 5k ) 

The matrix can be represented as 



V 



A(k) 



( 


1 





\ 




T 2 k 


T2k 










1 







T~2kTik 







\ 







T 5k J 



( rxk \ 

r 3fc r 3fc 

110 

10 1/ 



Let us examine the matrix 

/ rife 
r 2,fc+ir3fc 



C(k)B(k + 1) 



B{k)C{k). 
\ 



V 



r2,fc+ir 3 A : 
r2,fc+ir4 i fc+i r 2i fc+iT4 i fc + i T4 i fc + i 

<D 1 r 5ifc+ i r 5jfc+1 / 
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Taking into account that the matrix is triangular, we apply Theorem [5] 
to conclude that 

A = tr(E[C(l)B(l)]) = max(Eru, Er 21 + Er 3 i, Er 4 i, Er 51 ). 
6.4 A system with round routing 

Consider an open system depicted in Figure which consists of n = I + 1 
queues labeled with 0,1, ... ,1. Queue is intended to represent an external 
arrival stream of customers. Each incoming customer has to go to one of 
the other queues, being chosen by a regular round routing mechanism, and 
then leaves the system. 

With the routing mechanism, the customer that is the first to depart 
queue goes to the 1st queue, the second customer does to the 2nd queue, 
and so on. After the Ith customer directed to queue I, the next (I + l)st 
customer is directed to the 1st queue once again, and the procedure is 
further repeated round and round. 



1 




Figure 5: A system with round robin routing. 

The above system can be replaced with an equivalent fork-join network 
|21j which consists of n = 21 nodes (see Figure ED provided that for every 
nodes i = I + 1, 1 + 2, . . . , 21 in the new system, the service time is defined 
as 

T~ik = TQ^l-21+i- 

Suppose that 1 = 2. Evaluation of the system matrix gives 
A(k) = (I®T k GlfT k {I@G T l ) 



/ Tlk T lk T 3k Ti k T 3k \ 

T 2k T 2 fcT3fcT 4fc T 2 fcT3 fe T 4fc 

T 3k T 3k 

\ T 3k T ik T 3k T ik J 
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Figure 6: An equivalent fork-join network. 



Let us represent A(x) in the form 
A(k) 



( T lk T lk \ 

<D T~2k T 2 kT4k 

(DO 1 

\ T 4fc / 




B{k)C{k). 



Since the matrix product 
C{k)B{k + \) = 




<D r3fcr 4i fc+i 



has a triangular form, we get 

A = tr(E[C(l)B(2)]) = max(Er n , Er 2 i, Er 3 i + Er 4 i). 
Considering that ET31 = Er 4 i = Eroi , we finally have 
A = max(2Er i, Em, Et 2 i). 



7 Conclusion 

A stochastic dynamical system governed by the vector equation which is 
linear in some idempotent semiring was considered. New general conditions 
for Lyapunov exponent to exist for the system were given which involve 
the spectral radius of the mean state transition matrix. New method of 
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evaluation of the exponent was proposed based on a decomposition of the 
system state transition matrix. 

The above method was applied to the analysis of a class of queueing 
systems including open and closed tandem queues with finite and infinite 
buffers, fork-join networks, and systems with round-robin routing. Examples 
of evaluation of the mean service cycle time for the systems were considered, 
and related results in the form of some functions of the mean values of 
random variables that determine the service time were presented. 
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